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Abstract 

The process e"'"e~ W"'"W~7 is studied using the data collected by the L3 de- 
tector at LEP. New results, corresponding to an integrated luminosity of 427.4 pb~^ 
at centre-of-mass energies from 192 GeV to 207 GeV, are presented. 

The W"'"W~7 cross sections are measured to be in agreement with Standard 
Model expectations. No hints of anomalous quartic gauge boson couplings are 
observed. Limits at 95% confidence level are derived using also the process e"'"e~ — >■ 
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Introduction 



The increase of the LEP centre-of-mass energy well above the W boson pair-production thresh- 
old opens the possibility of studying the triple boson production process e"'"e~ — > W"'"W~7. We 
report on the cross section measurement of this inclusive process where the photon lies inside 
a defined phase space region. 

The three boson final state gives access to quartic gauge boson couplings represented by 
four-boson interaction vertices as shown in Figure^. At the LEP centre-of-mass energies the 
contribution of four-boson vertex diagrams, predicted by the Standard Model of electroweak 
interactions 0121, sue negligible with respect to the other competing diagrams, mainly initial- 
state radiation. The study of the W"*'W~7 process is thus sensitive to anomalous quartic gauge 
couplings (AQGC) in both the W^W~Z7 and W"''W~77 vertices. The presence of AQGC would 
increase the cross section and modify the photon energy spectrum of the W"'"W~7 process. This 
search is performed within the theoretical framework of References El and HI 

The existence of AQGC would also affect the e+e~ — >• uu^yj process via the W+W~ fusion 
diagram, shown in Figure^, containing the W^W~77 vertex The reaction e"'"e~ —>■ z/z/77 is 
dominated by initial-state radiation whereas the quartic Standard Model contribution from the 
W"'"W~ fusion is negligible at LEP. Also in this case the presence of AQGC would enhance the 
production rate, especially for the hard tail of the photon energy distribution and for photons 
produced at large angles with respect to the beam direction. 

The results are based on the high energy data sample collected with the L3 detector [Hj. 
Data at the centre-of-mass energy of y/s = 189 GeV, corresponding to an integrated luminosity 
of 176.8 pb~^, were already analysed and are used in the AQGC analysis. In the following, 
particular emphasis is given to the additional luminosity of 427.4 pb~^ recorded at centre-of- 
mass energies ranging from y/s = 192 GeV up to 207 GeV. 

The results derived on AQGC from the W"'"W~7 and uu'yy channels are eventually com- 
bined. 

Studies of triple gauge bosons production and AQGC were recently reported for both the 
W+W-7 ^ and Z77 H final states. 

Monte Carlo Simulation 

The dominant contributions to the W'^W~7 final state come from the radiative graphs with 
photons emitted by the incoming particles (ISR), by the decay products of the W bosons (FSR) 
or by the W's themselves (WSR). 

In this study, the signal is defined as the phase space region of the e"^e" W"'"W~7 process 
where the photon fulfills the following criteria: 

• > 5 GeV, where is the energy of the photon, 

• 20° < 6^ < 160°, where 6^ is the angle between the photon and the beam axis, 

• > 20°, where is the angle between the direction of the photon and that of the 
closest charged fermion. 

These requirements, used to enhance the effect of possible AQGC, largely contribute to 
avoid infrared and coUinear singularities in the calculation of the signal cross section. 
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In order to study efficiencies, background contaminations and AQGC effects, several Monte 
Carlo programs are used. 

The KORALW ^Oj generator, which does not include either the quartic coupling diagrams or 
the WSR, performs initial state multi-photon radiation in the full photon phase space. FSR 
from charged leptons in the event up to double bremsstrahlung is included using the PHOTOS [TT] 
package. The JETSET [12^ Monte Carlo program, which includes photons in the parton shower, 
is used to model the fragmentation and hadronization process. The KORALW program is used in 
the analysis for the determination of efficiencies. PYTHIA ^3] is used to simulate the background 
processes: e"'"e~ — >• Z/7 — >• qq(7), e"'"e~ — > ZZ — >• 4f(7) and e"'"e~ — >• Zee — >• ffee(7). 

The EEWWG ^ program is used to simulate the effect of AQGC. It includes 0{a) calculations 
for visible photons but is lacking the simulation of photons collinear to the beam pipe and of 
FSR. The net effect of collinear photons, included by implementing the EXCALIBUR jHj collinear 
radiator function, is to move the effective centre-of-mass energy towards lower values, reducing 
the expected signal cross section by about 18%. 

Other Monte Carlo programs which include WSR or full 0{a) corrections, such as YFSWW3 
[T3] and RacoonWW are used to cross check the calculations. 

For the simulation of the e^e~ uu'-f'-f process in the framework of the Standard Model 
the KORALZ jT2j Monte Carlo generator is used. NUNUGPV jlH] is also used to cross check the 
results, and found to be in agreement with KORALZ. The effects of AQGC are simulated using 
the EENUNUGGANO program . The missing higher order corrections due to ISR in EENUNUGGANO 
are also estimated by implementing the EXCALIBUR collinear radiator function. 

The response of the L3 detector is modelled with the GEANT ^Hj detector simulation 
program which includes the effects of energy loss, multiple scattering and showering in the 
detector materials and in the beam pipe. Time dependent detector inefficiencies are taken into 
account in the simulation. 

W+W~7 Event Selection and Cross Section 

The selection of W"''W~7 events follows two steps: ffist semileptonic W'^W^ qqeu or qq/xi^, 
and fully hadronic W+W~ qqqq events are selected [20] , then a search for isolated photons 
is performed. 

The photon identification in W"'"W~ events is optimized for each four-fermion final state. 
Photons are identified from energy clusters in the electromagnetic calorimeter not associated 
with any track in the central detector and with low activity in the nearby region of the hadron 
calorimeter. The profile of the shower must be consistent with that of an electromagnetic 
particle. Experimental cuts on photon energy and angles are applied, reflecting the phase 
space definition of the signal. 

Figure El shows the distributions of E^, 9^, and for the full data set, including the 
data at y/s = 189 GeV. Here, is defined as the angle of the photon with respect to the 
closest identified lepton or hadronic jet. Good agreement between data and Standard Model 
expectations is observed. Figure El shows the distributions of for the data collected at 

= 192 - 202 GeV and = 205 - 207 GeV, respectively. 

Table Q summarizes the selection yield. In total 86 W"'"W~7 candidate events are selected 
at y/s = 192 — 207 GeV. The Standard Model expectation, inside the specified phase space 
region, is of 87.8 ± 0.8 events. 

The quantity eww; representing the selection efficiencies for the W+W" — > qqlz/ and qqqq 
decay modes, ranges from 70% to 87%. The quantity is the photon identification efficiency 
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inside the selected phase space region. This efficiency takes into account small effects of events 
migrating from outside the signal region into the selected sample due to the finite detector 
resolution. Its value ranges from 52% to 80%, the lowest efficiencies being obtained in the fully 
hadronic sample where the high multiplicity makes the photon identification more difficult. 
The overall selection efficiency ^ww ^ ^7 is around 45% for all final states. 

The W"'"W~7 cross sections are evaluated channel by channel and then combined according 
to the Standard Model W boson branching fractions. The data samples at = 192 — 196 GeV, 
y/s = 200 — 202 GeV and y/s = 205 — 207 GeV are respectively merged. They correspond to the 
luminosity averaged centre-of-mass energies and to the integrated luminosities listed in Tabled 

The results, including the published value at ^/s = 189 GeV ^Tj, are: 

cTww7(188.6 GeV) = 0.29 ± 0.08 ± 0.02 pb (agM = 0.233 ± 0.012 pb) 

cTww7(194.4 GeV) = 0.23 ± 0.10 ± 0.02 pb (ctsm = 0.268 ± 0.013 pb) 

aww7(200.2 GeV) = 0.39 ± 0.12 ± 0.02 pb (ctsm = 0.305 ± 0.015 pb) 

cTww7(206.3 GeV) = 0.33 ± 0.09 ± 0.02 pb (ctsm = 0.323 ± 0.016 pb), 

where the ffist uncertainty is statistical and the second systematic. The measurements are 
in good agreement with the Standard Model expectations, (Tsm, calculated using EEWWG and 
reported with a theoretical uncertainty of 5% [21] ■ Figure |3] shows these results together with 
the predicted total W"^W~7 cross sections as a function of the centre-of-mass energy. 

The ratio between the measured cross section o"meas and the theoretical expectations is 
derived at each centre-of-mass energy. These values are then combined as: 

R = = 1.09 ± 0.17 ± 0.09 , 

O'SM 

where the ffist uncertainty is statistical and the second systematic. 

The systematic uncertainties arising in the inclusive W-pair event selections 120] are propa- 
gated to the final measurement and correspond to an uncertainty of 0.008 pb for all the energy 
points. Additional systematic uncertainties due to the electromagnetic calorimeter resolution 
and energy scale are found to be negligible. The total systematic uncertainty is dominated 
by the JETSET modelling of photons from meson decays (tt",?]). Its effect has been directly 
studied on data [7^ comparing the photon rate in e"'"e~ —>■ Z —>■ qq(7) events with Monte Carlo 
simulations. A correction factor of 1.2 ± 0.1 is applied to the rate of photons in the Monte 
Carlo simulation and its uncertainty is propagated. This uncertainty, fully correlated among 
the data taking periods, amounts to 6 % of the measured cross section. 



Determination of Anomalous Quartic Gauge Couplings 

The e+e" Process 

In the framework of References El and 01 the Standard Model Lagrangian of electroweak inter- 
actions is extended to include dimension-6 operators proportional to the three AQGC: a^/ h?, 
Oc/A^ and a^j ^ where A^ represents the energy scale for new physics. 

The two parameters a^j IS?' and Oc/A^, which are separately C and P conserving, generate 
anomalous W^W~77 and TIL'^^ vertices. The term an/h?, which is CP violating, gives rise 
to an anomalous contribution to the W^W~Z7 vertex. Indirect limits on ao/A^ and Oc/A^ 
were derived ^21; but only the study of W"'"W~7 events allows for a direct measurement of the 
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anomalous coupling anj 1^ ■ These couplings would manifest themselves by modifying the energy 
spectrum of the photons and the total cross section as shown in Figures El and IH respectively. 
The effect increases with increasing centre-of-mass energy. These predictions are obtained by 
reweighting the KORALW Monte Carlo events by the ratio of the differential distributions as 
calculated by the EEWWG and KORALW programs [7j. 

The derivation of AQGC is performed by fitting both the shape and the normalization of 
the photon energy spectrum in the range from 5 GeV to 35 GeV. Each of the AQGC is varied 
in turn fixing the other two to zero. 

The combination of all data, including the results at = 189 GeV [7 , gives: 

ao/A^ = 0.000 ± 0.010 GeV^^ 
ajk^ = -0.013 ± 0.023 GeV^^ 
a„/A2 = -0.002 ± 0.076 GeV^^ 

where systematic uncertainties are included. 

At the 95% confidence level, the AQGC are constrained to: 

-0.017 GeV"2 < ao/A^ < 0.017 GeV~2 
-0.052 GeV-^ < aJA"^ < 0.026 GeV"^ 
-0.14 GeV~2 < an/A^ < 0.13 GeV^^ 

All these results are in agreement with the Standard Model expectation. The sign of the Oq 
and AQGC, obtained with the EEWWG reweighting, is reversed according to the discussions 
in References I2SI and 1211 

The e+e" — > l/z>77 Process 

The sensitivity of the e"'"e~ — > z/z/77 process to the ao/A^ and ad A^ AQGC, through the dia- 
grams shown in Figure^, is also exploited. Events with an acoplanar multi-photon signature 
are selected [23]. In this letter we report on results from data at ^/s = 192 — 207 GeV. 

Figure El shows the two-photon recoil mass, Mrec, distribution, with the predicted AQGC 
signal for a non-zero anomalous coupling oq/A^. The number of selected events in the Z peak 
region, defined as 75 GeV < Mrec < 110 GeV, is 43 in agreement with the Standard Model 
expectation of 47.6 ± 0.7. 

The AQGC signal prediction is reliable only for recoil masses lower than the mass of the Z 
boson as the interference with the Standard Model processes is not included in the calculation. 
Requiring Mrec < 75 GeV, no event is retained by the selection in agreement with the Standard 
Model expectation of 0.35 ± 0.05 events. 

A reweighting technique based on the full matrix elements as calculated by the EENUNUGGANO 
Monte Carlo, is used to derive the AQGC values. Including the results at ^/s = 183 GeV and 
^/s = 189 GeV [7j, the 95% confidence level upper limits are: 

-0.031 GeV"^ < ao/A^ < 0.031 GeV"^ 
-0.090 GeV"^ < aJA^ < 0.090 GeV^^ , 

where only one parameter is varied at a time. 

The dominant systematic uncertainty comes from the theoretical uncertainty of 5% [21] in 
the calculation of anomalous cross sections. 
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Combined Results 

The results obtained from the W+W~7 and uu^yj processes are combined. No evidence of 
AQGC is found and 95% confidence level limits are obtained separately on each coupling as: 

-0.015 GcV^2 < ao/A^ < 0.015 GeV^^ 
-0.048 GeV-2 < aJA^ < 0.026 GeV"^ 
-0.14 GeV-2 < a„/A2 < 0.13 GeV'^. 

Appendix 

The results on the e+e^ W+W^7 cross sections are also expressed for a different phase space 
region defined by: 

m E^> 5 GeV, 

• |cos(^^)| < 0.95, 

• cos(a„J < 0.90, 

• Mffi = Mw ± 2rw, where Mff are the two fermion-pair invariant masses. 
The results read: 

C7WW7(188.6 GcV) = 0.20 ± 0.09 ± 0.01 pb (ctsm = 0.190 ± 0.010 pb) 

aww7( 194.4 GeV) = 0.17 ± 0.10 ± 0.01 pb (ctsm = 0.219 ± 0.011 pb) 

cTww7(200.2 GeV) = 0.43 ± 0.13 ± 0.02 pb (ctsm = 0.242 ± 0.012 pb) 

aww7 (206.3 GeV) = 0.13 ± 0.08 ± 0.01 pb (cjsm = 0.259 ± 0.013 pb), 

where the first uncertainty is statistical, the second systematic and the values in parentheses 
indicate the Standard Model predictions. 
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' data 






ATexp 
1^ TOT 




194.4 GeV 


qqei/7 


4 


0.748 ± 0.007 


0.613 ± 0.023 


3.76 ± 0.14 


1.41 ± 0.09 


(113.4 pb-i) 


qq/iZ/7 


6 


0.731 ±0.007 


0.728 ±0.026 


4.52 ±0.15 


1.68 ±0.09 




qqqq7 


9 


0.865 ± 0.004 


0.537 ±0.013 


12.89 ± 0.34 


5.22 ±0.28 


200.2 GeV 


qqei/7 


5 


0.726 ± 0.007 


0.631 ± 0.023 


4.28 ±0.14 


1.47 ±0.08 


(119.8 pb-i) 


qq//i/7 


4 


0.712 ± 0.007 


0.744 ± 0.025 


5.50 ±0.16 


2.07 ±0.10 




qqqq7 


18 


0.836 ± 0.004 


0.521 ±0.012 


13.71 ± 0.29 


5.31 ±0.28 


206.3 GeV 


qqez/7 


7 


0.700 ± 0.005 


0.626 ± 0.024 


6.58 ±0.22 


2.21 ±0.13 


(194.2 pb-^) 


qq//i/7 


4 


0.714 ±0.005 


0.787 ±0.026 


9.54 ±0.26 


3.48 ±0.16 




qqqq7 


29 


0.823 ± 0.005 


0.540 ±0.011 


26.97 ±0.50 


12.11 ±0.39 



Table 1: Number of observed events, N^atai W+W~ and W+W~7 selection efficiencies, £ww and 
expected total number of events, N^qj., and background estimates, N^Bkgri f*^^ various 
decay channels according to the Standard Model prediction. The background estimates include 
FSR and misidentified photons. All uncertainties come from Monte Carlo statistics. The 
average centre-of-mass energies, and the integrated luminosity of the three subsamples 

are also listed. 
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Figure 1: Feynman diagrams containing four-boson vertices leading to the a) W^W 7 and b) 
i/p77 final states. 
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Figure 2: Distributions at the centre-of-mass energies a/s = 189 — 207 GeV of a) the photon 
energy, b) the cosine of the angle of the photon to the beam axis and c) the cosine of the angle 
of the photon to the closest charged lepton or jet. 
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Figure 3: Distributions of the photon energy for the semileptonic qqev, qqjj.v and fully hadronic 
W+W~7 decay modes corresponding to the data collected at a) ^/s — 192 — 202 GeV and b) 
^ — 205 — 207 GeV. The cross-hatched area is the background component from WW, ZZ, 
Zee, and qq(7) events. The FSR distribution includes the contribution of photons radiated 
off the charged fermions and photons originating from isolated meson decays. Distributions 
corresponding to non-zero values of the anomalous coupling a„/A^ are shown as dashed lines. 
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Figure 4: Measured cross section for the process e"'"e~ ^^^^^^"7 compared to the Standard 
Model cross section as a function of the centre-of-mass energy, as predicted by the EEWWG 
Monte Carlo within phase-space requirements. The shaded band corresponds to a theoretical 
uncertainty of ±5%. The three dash-dotted lines correspond to the cross section for the 
indicated values of the anomalous coupling a„/A^ (in GeV~^ units). 
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Figure 5: Recoil mass spectrum of the acoplanar photon pair events selected at ^/s — 192 — 
207 GeV. 
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